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A NEAR-INFRARED STUDY OF THE HIGHLY-OBSCURED ACTIVE 

STAR-FORMING REGION W51B 
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ABSTRACT 

We present wide-field JiJA s -band photometric observations of the three compact H II regions G48.9- 
0.3, G49.0-0.3, and G49.2-0.3 in the active star-forming region W51B. The star clusters inside the three 
compact H II regions show the excess number of stars in the J—K s histograms compared with reference 
fields. While the mean color excess ratio (Ej_h / Eh-kJ of the three compact H II regions are similar 
to ~ 2.07, the visual extinctions toward them are somewhat different: ~ 17 mag for G48.9-0.3 and 
G49.0-0.3; ~ 23 mag for G49.2-0.3. Based on their sizes and brightnesses, we suggest that the age of 
each compact H II region is < 2 Myr. The inferred total stellar mass, ~ 1.4 x 10 4 M©, of W51B makes 
it one of the most active star forming regions in the Galaxy with the star formation efficiency of ~ 10 



%. 
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I. INTRODUCTION 

W51B is a giant molecular cloud (GMC) complex 
located at the tangential point of the Sagittarius spi- 
ral arm of the Galaxy, and it is a very active star- 
forming region that consists of several compact H II re- 
gions. In radio continuum and molecular line images, 
the compact H II regions appear to be still embedded in 
parental molecular clouds, indicating that they are rela- 
tively young (e.g. Bieging 1975; Koo & Moon 1997; Koo 
1997, 1999; Moon & Park 1998). The H II regions also 
show strong X-ray emission likely caused by magnetic 
activities of internal young stars (Koo et al. 2002). The 
line-of-sight velocities (61-70 km s^ 1 ) of the sources in 
W51B are greater than the maximum kinematic ve- 
locity (54-58 km s _1 ) permitted by the Galactic ro- 
tation (Koo 1999, and references therein). This has 
been known as the "high-velocity stream" in this re- 
gion which might have resulted from the spiral density 
wave associated with the Sagittarius spiral arm. The 
southern part of W51B appears to show an expanding 
ring structure with two molecular clumps associated 
with the compact H II regions G48.9-0.3 and G49.0- 
0.3 (Moon & Park 1998). Using 21-cm H I absorption 
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spectra toward the compact H II regions, Kolpak et al. 
(2003) determined the distance to W51B to be ~ 5.6 
kpc.* 

One of the most interesting features of W51B is its 
very high star formation efficiency (SFE): while the 
Galactic average value is 1-2 % (Myers et al. 1986), 
Moon & Koo (1994) obtained 16 % for the overall SFE 
of W51B and Koo (1999) did 7 % and 15 % for the 
north-eastern and south-western parts of W51B, re- 
spectively. Most of the previous studies of the star 
formation in W51B, including the ones carried out by 
Moon & Koo (1994) and Koo (1999), have exclusively 
used radio contiuum and molecular line data. Recently, 
Nanda Kumar et al. (2004) have obtained several small- 
field (= 90") near-infrared (NIR) images of W51B, and 
identified a new, small-scale clumpy structure of star 
clusters of the compact H II regions. However, the 
observations of Nanda Kumar et al. (2004) were some- 
what shallow, only reaching ~16.2, ~14.4, and ~15.8 
mag for J-, H-, and A-band, respectively. This rather 
shallow exposure, together with their small field of 
view, made it difficult to study fainter, extended stel- 
lar populations in W51B. They calculated the mass of 
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*We will use 5.6 kpc as the distance to W51B in this paper. 
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Right Ascension ( J2000) 



Fig. 1. — Gray-scale map of the CO J=l-0 line emission 
integrated over « LSR between +45 and +75 km s _1 (Car- 
penter & Sanders 1998), superposed on the black contours 
of the 21-cm radio continuum map (Koo & Moon 1997). 
The contour levels are 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 
0.9, and 1.0 Jy beam -1 . The white boxes and circles rep- 
resent the NIR observation fields and the NIR (especially 
-Ks-band) nebulous area (see Figure 5): SW-Field includ- 
ing G48.9-0.3 and G49.0-0.3 (from SW to NE), and NE- 
Field including G49.2-0.3. The peak positions of the radio 
continuum are consistent with those of the NIR nebulosity. 

each cluster assuming an average visual extinction of 7 
mag. With such moderate extinction, however, all of 
OB-type members should be visible in second Palomar 
Observatory Sky Survey (POSS II) red image (i? Cj ii m 
= 20.8 mag, Reid et al. 1991), while no optical coun- 
terpart is known. This implies therefore the extinc- 
tion toward W51B H II regions must be larger than 7 
mag. In this paper, we obtained deeper NIR photo- 
metric data that cover much more extended regions of 
W51B. With these new data, we find the initial mass 
function and the total stellar mass of star clusters in 
W51B, as well as NIR-based SFE of W51B. Our new 
observations confirm the results of the previous radio 
observations that W51B is one of the most active star 
forming regions in the Galaxy. 

II. OBSERVATIONS 

We carried out near-IR J H K s -b&nd photometric ob- 
servations of W51B using Wide-field InfraRed Cam- 
era (WIRC) aboard the 5-m Palomar Hale telescope 
on 2003 June 17 and 18. WIRC is equipped with a 
HAWAII-II HgCdTe 2K array and the pixel scale is 
0.25" pixel -1 , providing a 8.7'x 8.7' field of view. The 



typical seeing was 1" over the observations. We selected 
two fields based on the previous radio observations, and 
Figure 1 shows the locations of the observed fields su- 
perimposed on previous radio images. As in Figure 1, 
one field includes G48.9-0.3 and G49.0-0.3, while the 
other field contains G49.2-0.3, the brightest radio con- 
tinuum source in W51B (Koo & Moon 1997), at its 
center. We shall call the former SW-Field and the lat- 
ter NE-Field hereafter. We obtained 30 and 60 dithered 
frames of SW-Field and NE-Field respectively, for each 
JHK S band. The exposure of individual frame was 10 
sec. Each dithered frame was shifted and combined af- 
ter subtracting dark and median-combined sky frames 
and correcting flat field. Figures 2 and 3 present the 
final three-color images of the two fields, where we can 
identify extended emission around the H II regions due 
to Br-7 emission in the K s band. 

In order to find stars in our WIRC images, we itera- 
tively used the PSF-fitting program in IRAF t /DAOP- 
HOT until the residuals show no apparent stars left in 
the images. As a result, we have detected > 11000 and 
> 7000 stars in the SW-Field and NE-Field, respec- 
tively. Among them ~ 10000 and ~ 6000 stars have 
signal-to-noise ratio greater than 10 in all the JHK S 
bands, within 10-er limiting magnitudes of J = 19.8, 
H = 19.0, and K s = 17.9 mag. For photometric and 
astrometric solutions, we used ~ 200 stars of 2MASS 
Point Source Catalog* in our fields. The magnitudes 
of the faintest stars were J = 21.5 ± 0.5, H = 19.7 ± 
0.2, and K s = 19.2 ± 0.4 for SW-Field and J = 21.5 
± 0.5, H = 19.8 ± 0.2, and K s = 19.5 ± 0.5 for NE- 
Field respectively, while the magnitude of the brightest 
and unsaturated star was ~ 10 in all bands. Our as- 
trometric solutions gave 1-a uncertainties of ~0"3 for 
SW-Field and ~07l for NE-Field with respect to the 
2MASS system. 

III. Star Clusters in W51B 
(a) Membership Selection 

In order to investigate the existence of noticeable 
star clusters inside the H II regions, we first examine 
the stellar number density map as in Figure 4 which 
was obtained by Gaussian smoothing of number of stars 
with 15" FWHM. It shows clusterings of stars inside 
the three compact H II regions. Although there are 
many stars at south-east of NE-Field, they should not 
be regarded as cluster members because they are ran- 
domly distributed outside the cloud. They may be 
background stars behind W51B GMCs. In order to de- 
termine the cluster members and examine their proper- 
ties, we define the cluster regions and the corresponding 



^IRAF (Imaging Reduction & Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science 
Foundation. 

t http : / / www. ipac.caltech.edu/2mass/ 
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Fig. 2. — J, H, and K a three-color image of SW-Field, 
including the two compact H II regions G48.9-0.3 and 
G49. 0-0.3. The showing field is 8.2' x 8.2' centered at 
(aj2ooo.o, <5j2ooo.o) = (19 h 22 m 21 s , +14°4'26"). The J-, H-, 
and A's-band data are represented as blue, green, and red, 
respectively. 




Fig. 3. — Same as Figure 2 but for NE-Field, which in- 
cludes G49.2-0.3. The showing field is 6.7' x 7.0' centered 
at O./2000.0, <5j2ooo.o) = (19 h 23 m 7 s , +14°16'47"). 



reference regions on the basis of the positions of H II 
regions, which are described as solid and dashed cir- 
cles, respectively, in Figure 5. Each reference region is 
selected just outside the cluster region to minimize the 
offset effect of spatial density fluctuation. A reference 
region has the same area with the corresponding clus- 
ter region. There are about 600, 1500, and 500 stars in 
the cluster regions within the H II regions G48.9-0.3, 
G49.0-0.3, and G49.2-0.3, respectively. These numbers 
of stars inside the cluster regions are at least 1-er larger 
than the counts in the randomly selected regions hav- 
ing same area. For G49.2-0.3, we avoid selecting the 
random regions outside molecular cloud because there 
must be a lot of background field stars in the off-cloud 
area. 

Figure 6 shows the histograms of J—K s color of the 
stars detected at all bands in the direction of three H II 
regions. The first bump at J—K s ~ 1.4 is found on the 
histograms for both cluster and reference regions, thus 
it should be caused by foreground field stars. With an 
assumption of Ay ~ 1 mag kpc , they seem to be 
mainly composed of the Sagittarius arm which extends 
through a range of ~ 3-6 kpc at / ~ 49° (Wainscoat et 
al. 1992). We see outstanding bumps toward G48.9-0.3 
and G49.0-0.3 at 2.2 < J-K s < 3.2, compared to the 
histograms of the reference regions. In the histogram 
for G49.2-0.3, the extra bump is at 2.6 < J-K s < 4.5. 
The inset panel represents the excess number of stars 
compared to reference region. The excess should be 
caused by the probable cluster members. 

Although the bump at J—K s > 2.2 probably in- 
cludes most of the cluster members (we call them 
"B" (Bump) members hereafter), we may fail to se- 
lect the pre-main-sequence (PMS) stars which are gen- 
erally redder than the main-sequence (MS) stars due 
to the intrinsic emission from the circumstellar dust. 
Thus, we need to separately include the PMS stars 
into the cluster member candidates. We use H—K s vs. 
J—H color-color diagram (CCD) to select PMS stars 
with intrinsic NIR color excess. In order to distinguish 
the intrinsic color excess from the reddening by inter- 
stellar extinction, we first compute the color excess ra- 
tio Ej- H /Eh-k s toward W51B in the 2MASS system 
by a least-squares fit with 3-cr clipping. Figure 7 shows 
the CCD of bright stars (Kg < 13) with each color error 
less than 0.01 mag. The stars after 3-cr clipping yield 
the ratios of 2.08 ± 0.04 for SW-Field and 2.07 ± 0.02 
for NE-Field. Those are similar but slightly larger than 
1.98 of Fitzpatrick (1999) provided from extensive data 
in the Galaxy or 1.97±0.01 of Han (2001) obtained for 
a 25'x25' field including whole W51B giant molecular 
cloud region and large off-cloud region. Because our 
fields are focused on the cloud core regions, there could 
be a little difference from the value calculated for a wide 
field. Hereafter, we apply the color excess ratio of 2.07 
in the direction of the H II regions in W51B and use 
the relations of Fitzpatrick in order to estimate total 
extinction. 
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Fig. 4. — Gray-scale map of the number density of stars detected in all bands, superposed on red and blue contours of 
the number densities of "P" and "B" members, respectively (see text), for each field. The colored thick and thin contours 
represent 3-cr and 2-cr of the number density where a is the standard deviation of each kind of density. The white contours 
are the 21-cm radio continuum whose levels are given in Figure 1. The left and right panels are for SW-Field and NE-Field, 
respectively. 
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Fig. 5. — Tfs-band images for SW-Field (left) and NE-Field (right). For NE-Field, the observed J-, H-, and K B -band 
frames were shifted each other, and we use only the matched area which is inside the black box. The solid and dashed circles 
represent the cluster and reference regions, respectively. The reference region has the same area with the corresponding 
cluster region. Also shown are the cluster-core regions which are defined inside the 3-cr stellar density contours in Figure 4. 
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Fig. 6. — Histograms along J—K s color. Solid and dot- 
ted lines show the histograms for a cluster region and the 
corresponding reference region, respectively. The error bar 
denotes the Poisson noise level in each bin. The difference 
for the second bump in the number of stars between the 
cluster region and the reference region is displayed as an in- 
set of each panel. The white and black histograms represent 
the numbers of PMS and the others, respectively. 
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Fig. 7. — Color-color diagram of the bright stars {K B <12>) 
in each field with good photometry (H—K s and J—H color 
errors less than 0.01 mag). Reddening slopes are obtained 
using a least-squares fit with 3-ct clipping. The left panel is 
for SW-Field, and the right one for NE-Field. The line in 
each panel shows the fitted reddening slope, both of which 
are ~ 2.07. 



Figure 8 is H—K s vs. J—H CCDs where the red- 
dening slope of the dashed lines is determined by the 
color excess ratio obtained in the previous paragraph. 
The area between the two dashed lines is the "redden- 
ing zone" in which the reddened dwarfs and giants fall. 
Stars with intrinsic infrared excess emission can fall into 
the right of the reddening zone. Lada & Adams (1992) 
called this region as "excess zone" . The sources in the 
excess zone are generally pre- and post-main-sequence 
stars with circumstellar dust and galaxies with redshifts 
z > 0.1. Toward the compact H II regions, most of them 
should be pre-main-sequence stars. We mark these 
stars in red color, except possible spreading foreground 
stars with J-K s < 2.2. We call them "P" (PMS) mem- 
bers hereafter. And we mark the "B" members in blue 
color. Most of the "B" members are MS stars, though 
some of them are marked in both colors which indicate 
that they are probably PMS stars. Since most stars in 
these categories are located above the reddening line 
of 3M stars in J—K s vs. J color-magnitude diagrams 
in Figure 9, we suppose the observed cluster members 
to be OB-type MS stars and/or intermediate- to high- 
mass PMS stars. We note the numbers of stars in the 
two categories in Table 1. The selected cluster member 
candidates are 168, 589, and 120 for G48. 9-0.3, G49.0- 
0.3, and G49.2-0.3, respectively. These candidates, 
however, include some field star contamination though 
we already excluded most of the foreground field stars 
composing the bump at J—K s < 2.2 in Figure 6. 

In order to analyze how many field stars are con- 
tained, we choose out the stars inside the reference 
region in each category using the same method, and 
compare the number of stars between cluster and ref- 
erence regions. We confirm that at least 60-80 % of the 
selected member candidates are statistically reliable for 
G48.9-0.3 and G49.2-0.3, while the reliability is 40 % 
for G49.0-0.3. Since G49.0-0.3 is sparsely distributed 
and the extinction varies a lot, the cluster members 
could not be determined more reliably. However, most 
of the field star contamination contributes to the low lu- 
minosity regime (J > 18 mag), thus their contribution 
to the total luminosity is expected to be insignificant. 

(b) Interstellar Extinction toward W51B Clus- 
ters 

We expect that the average interstellar extinction 
toward the cluster can be simply obtained from the 
mean photometric color of members. Since the range 
of the intrinsic NIR colors of OB-type MS stars is very 
narrow (A(J — K s ) intrinsic ~ 0.2 mag) compared to 
the J—K s range of sources in the observed field, the 
J—K s color can be a good indicator of extinction. 

We transform the observed J—K s colors of MS mem- 
ber candidates to the visual extinctions using 

A V = 6.31 [(J - A s ) obs . - (J - intrinsic], (1) 

which is obtained from Fitzpatrick's relations converted 
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Table 1 

Number of Cluster Member Candidates 



H II region "B" members "P" members "B" & "P" members 



G48. 9-0.3 143 66 41 

G49.0-0.3 560 105 76 

G49. 2-0.3 109 45 34 



Note. — The fourth column represents the number of stars in both 
categories, and we include these stars into PMS member candidates. 
MS member candidates are considered as the stars in "B" but not "P" 
category. 



to the 2MASS system^ which is based on Carpenter 
(2001). We assume the intrinsic J—K s color of OB 
stars to be approximately zero. PMS stars must pos- 
sess redder intrinsic colors though there is not a defi- 
nite model explaining the color of massive PMS stars 
yet. We assume the average extinction of PMS stars 
is identical to that of MS stars. The mean visual ex- 
tinctions for the member candidates are 17 mag for 
both G48. 9-0.3 and G49. 0-0.3, and 23 mag for G49.2- 
0.3. Similarly, mean J-band extinctions are obtained 
from Aj = 1.73 (J - K s ) ohs .: 4.6 mag for G48.9-0.3 
and G49.0-0.3, and 6.2 mag for G49.2-0.3. The photo- 
metric diagrams of G49.2-0.3 show a wide separation 
between the cluster and foreground stars at J—K s ~ 
2, which is consistent with the higher extinction value 
than the other two clusters. 

(c) Spatial Distribution of Member Stars 

In order to see the clustering feature of the member 
candidates rigorously, we describe the number densities 
of the "P" and "B" members as red and blue contours 
in Figure 4, overlapped on the gray-scale map for the 
number density of all observed stars. We define the 
cluster-core region inside the contour of 3-cr where a 
is the standard deviation of the density fluctuations. 
In Table 2, we display the stellar column and volume 
densities of each group inside each cluster-core region 
after statistical subtraction of the stars outside the re- 
gion under consideration. They are the lower limits 
of the cluster densities because the lowest mass of the 
marginally detectable members are 2-5M correspond- 
ing to the 10-<7 limiting magnitudes. 

In G48.9-0.3, both MS and PMS members are well 
confined inside the area of radio continuum, elongated 
from north to south. The PMS members appear to 
be peanut-shaped. G49.0-0.3 shows a large and bright 
density peak of "B" members. The PMS stars are con- 
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centrated on the southern part of it near the radio con- 
tinuum peak. There are two small PMS clusters at 
rather distant places (~ 2') on the north-west of G49.0- 
0.3. The central region of G49.2-0.3 is crowded inside a 
small place, while the radio continuum source is much 
larger than the previously mentioned two other H II re- 
gions. If G49.2-0.3 is a "blister-type" H II region being 
suggested from the radio continuum morphology and 
the distribution of molecular gas (e.g. Koo 1999), the 
large size of the ionized gas could be explained by low 
density of ambient medium. 

IV. STAR FORMATION EFFICIENCY 

With the J-band magnitudes of member stars, we 
find the initial mass function (IMF) and the total stellar 
mass of each cluster. We use J-band luminosity which 
is less affected by the intrinsic excess emission of PMS 
stars than H- or if s -band luminosity. We transform 
the observed J magnitudes of all member candidates 
to the absolute magnitudes applying the mean extinc- 
tion of the member candidates derived at § III (b). The 
transformation from the absolute J-magnitude to the 
initial mass is accomplished using the theoretical mass- 
luminosity relation from the Geneva models and theo- 
retical colors of MS stars (Schaller et al. 1992; Bessell 
et al. 1998) and Testi et al. (1998) for PMS stars less 
massive than 3M Q . Since massive PMS stars (> 3M Q ) 
would be located very close to the main sequence in the 
mass-luminosity relation, the MS relation is applied for 
them. We assume the age of clusters to be 1 Myr, and 
the derived total mass is not sensitive to the age. Figure 
10 shows the resultant IMF for each cluster. The num- 
ber of stars increases as the stellar mass decreases at 
the high mass regime, while the detection limit causes 
the IMFs to decrease at the lower mass regime. For 
comparison, we plot the IMF of the young open cluster 
NGC 2264 (Sung et al. 2004, solid curve) extrapolated 
using the IMF of the solar neighborhoods (Lee & Sung 
1995) for the massive stars (> 30M Q ) and vertically 
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Table 2 

Number of Stars in Each Cluster-Core Region 



H II region 


Name a 


Number 


Area (□") 


N/A pc2 b 


N/V pc s c 


G48.9-0.3 


P4855044-0017085 
B4855164-0016594 


36 ± 7 
23 ± 6 


2450 
449 


63 ± 12 
218 ± 55 


20 ± 4 
160 ± 40 


G49.0-0.3 


P4859241-0018016 
B4859462-0017400 


8 ± 4 
122 ± 16 


1125 
3674 


30 ± 16 
142 ± 18 


14 ± 7 
36 ± 5 


G49. 2-0.3 


P4912136-0020508 
B4912111-0021025 


16 ± 5 
26 ± 7 


1989 
1064 


34 ± 10 
105 ± 27 


12 ± 4 
50 ± 13 



a The position of the local peak in Galactic coordinates were used to name the clustering groups. 
The first letter represents "PMS" and "Bump" , respectively. 

b Column density. N/A pc 2 = 4262 N/A^n adopting the distance of 5.6 kpc. 
c Volume density. N/V pc 3 = 66430 N/A 3 '?, if we assume the cluster is spherical. 



shifted to fit our data within 1.0 < log (M/M Q ) < 1.3. 
Sung ct al. noticed that NGC 2264 has similar shape of 
IMF to the Pleiades and the Trapezium. The IMFs of 
G48.9-0.3 and G49.0-0.3 agree well with the compari- 
son IMF at M > 5-6 M . For G49.2-0.3, the consis- 
tent range is M > 10M Q due to high extinction toward 
this cluster. 

In order to estimate total mass of a cluster, the uni- 
versality of IMF is assumed. The total mass is deter- 
mined from 



M t , 



J M£ (log M) d log M (2) 



where £ is the IMF of NGC 2264 which described in the 
previous paragraph. T he upper limit of the integral 
is decided on the mass of the most massive stars, and 
the lower limit is fixed to O.OIMq. Table 3 represents 
the total mass of the cluster associated with each H II 
region. We also list the stellar mass of the only ob- 
served stars, instead assuming the universality of IMF. 
Undetected stars would contribute to the cluster mass 
in amount of at least twice than the observed massive 
stars do. 

We calculate the star formation efficiency (SFE) 
which is defined as the ratio of the stellar mass to 
the total mass of stars plus cloud, using total stellar 
mass in Table 3 and the molecular cloud mass in Ta- 
ble 1 of Koo (1999). The overall SFEs are 7 % for the 
cloud SW including G48.9-0.3 and G49.0-0.3, and 17 
% for the cloud NE including G49.2-0.3. Even if we 
only use the stellar mass of the observed massive stars, 
the efficiencies are estimated to be 3-4 %. Because 
no cluster is composed of only massive OB stars, the 
SFEs of W51B appear to be clearly high compared to 
the Galactic average value of 1-2 %. Besides, the high 



SFEs are lower limits because we ignore the possible 
smaller clusters/associations along the GMCs. For ex- 
planations about the high SFE, a cloud-cloud collision 
between high-velocity (68 km s _1 ) cloud and the W51 
GMC (Carpenter & Sanders 1998) and the enhance- 
ment due to the spiral density wave (Koo 1999; Nanda 
Kumar et al. 2004) have been speculated. 

Table 2 of Lada & Lada (2003) shows that the SFEs 
for nearby embedded clusters range from ~ 10 % to 
30 % when the cloud core masses are adopted, not the 
masses for entire giant molecular clouds. The com- 
pressed core masses presented in that paper are less 
than 1000 M , which is at least 10 times less than the 
entire mass of W51B GMCs. We cannot compare the 
core SFEs to the overall SFEs of W51B GMCs until the 
individual core masses are obtained by high resolution 
observations. 

V. AGE OF CLUSTERS 

As the stars emit Lyman continuum photons, the 
H II region increases in size. Therefore, in a given size 
of the H II region, the age is inversely proportional to 
the Lyman continuum photon flux. If the H II region 
expands in a uniform medium, the radius of an H II 
region as a function of time is given by 



4/7 



where 



3N r 



1/3 



4irn e riH/3 



(3) 



(4) 



is the initial Stromgren radius, a is the sound speed 
in the ionized gas (~ 10 km s _1 for 10 4 K), n e and 
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Table 3 

Total Stellar Mass Estimations 



H II rej 


,'ion 


Method 


Mass Range 


Total Stellar Mass 








(M 


) 


(Mo) 


G48.9- 


-0.3 


1 


0.01 ~ 


20 


2400 






2 


1.3 ~ 


20 


800 


G49.0- 


-0.3 


1 


0.01 ~ 


25 


4200 






2 


0.8 - 


25 


2300 


G49.2 


-0.3 


1 


0.01 ~ 


100 


7800 






2 


2.5 - 


100 


1300 



Note. — For the method (1), we assume the combined IMF of 
NGC 2264 and solar neighborhoods as a universal IMF in order to 
estimate total mass including undetected low-mass stars. Method 
(2) is the stellar mass of the observed stars without the assumption. 




Fig. 8. — Color-color diagram of stars in each cluster region. The member candidates are marked in red and blue colors. 
The red marks are "P" members, while the blue marks are the "B" members. The red open circles and the blue closed 
circles represent the member candidates with small color error (< 0.1 mag), and the colored dots are the candidates with 
larger error. The thin and thick solid curves are the loci occupied by the unreddened dwarfs and giants (Bessell & Brett 
1988; Koornneef 1983). The dashed lines are parallel to the reddening vector. 



NIR STUDY OF W51B STAR-FORMING REGION 



25 




Fig. 9. — Color-magnitude diagram of stars in each cluster region. The black, thick curve represents the unreddened 
zero-age main-sequence (ZAMS) relation for Z=0.020 from the Geneva models (Schaller et al. 1992; Bessell et al. 1998), 
adopting the distance of 5.6 kpc to the compact H II regions. The dashed lines are parallel to the reddening vector of ZAMS 
star having mass of 120M© (up) and 3Mq (down). The black, thin curves are reddened ZAMS corresponding with Ay of 5, 
10, 15, 20, 25, and 30 mag. Additionally, the red curve is the loci of unreddened PMS stars with age of lMyr (Testi et al. 
1998). The red and blue marks represent the member candidates which are described in Figure 8. 




Fig. 10. — Initial mass functions for all cluster member candidates (black) and the member candidates with color error 
less than 0.1 mag (red). We count the stars within the mass interval A log (M/Mq) — 0.2 and in the same-sized bin shifted 
by 0.1 to reduce the binning effect. They are represented as open and closed circles. The error bars denote the statistical 
fluctuations in the counts under the assumption that the Poisson statistics governs the fluctuations. The initial mass function 
of NGC 2264 (solid curve) is vertically shifted to fit our data and used to estimate total mass of W51B. 
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Table 4 

Properties of the Clusters Associated with Three H II regions 



H II region 


(Ay) 


Radius 


Spectral Type of 


\ogN c a 


Age b 




(mag) 


(PC) 


the most massive star 




(Myr) 


G48.9-0.3 


17 


1.6 


BO 


47.85 (46.90) 


0.6 (1.1) 


G49.0-0.3 


17 


2.2 


09 


48.31 (48.23) 


0.8 (0.9) 


G49. 2-0.3 


23 


2.7 


04 


50.43 (50.25) 


0.3 (0.3) 



a Lyman continuum photon flux of all members (and the members without NIR excess). 
b Dynamic age of the H II region using all members (and the members without NIR excess). 



nn are the ambient densities of electron and atomic 
hydrogen, (3 is the recombination coefficient (~ 3 x 

10" 13 cm 3 s- 1 for T e = 10 4 K), and N' c is the effective 
Lyman continuum photon flux (e.g. Dyson & Williams 

1980). N c is evaluated assuming 90 % of the Lyman 
photons are absorbed by dust within the H II region 
(Wood & Churchwell 1989). If we adopt n e — uh, the 
age of an H II region is simplified as follows: 




Although the three H II regions in W51B are not 
spherical, we may roughly estimate the age of the H II 
regions from equation (5). The mean electron density 
of W51B H II regions is ~ 10 3 cm -3 (Mufson & Liszt 
1979; Mehringer 1994). The present sizes of the H II 
regions are decided from the appearance of the 21-cm 
radio continuum emission (0.3 Jy beam -1 contour in 
Figure 1), and are presented in Table 4. In the table, 
we also present the spectral type of the most massive 
star of a cluster and the total Lyman continuum pho- 
ton flux that the cluster members emit at ZAMS stage 
(Panagia 1973), together with the resultant age. The 
Lyman flux and the inferred age using only the stars 
which are expected to be already in main-sequence are 
put in parentheses. We derive the age to be the order 
of or younger than 1 Myr for all H II regions. However, 
even in an ideal case, the dynamic age of the H II re- 
gion may give only the age of the most massive stars 
of the cluster. In addition, the dynamic age is quite 
sensitive to the size of H II region, which depends on 
the distance. Nevertheless, the deduced dynamic age is 
roughly consistent with the age estimation via JHK S 
excess fraction, which is shown in next paragraph. 

The fraction of the stars possessing NIR excess emis- 
sion is revealed to be another indicator of age (Haisch 
et al. 2001a, and therein). For a very young embedded 



cluster NGC 2024 (~ 0.3 Myr, Meyer 1996), the JHK 
excess fraction of 58% ± 7% is derived by Haisch ct 
al. (2000). The fraction decreases as the age increases: 
50% ± 7% (Lada et al. 2000) and 21% ± 5% (Haisch et 
al. 2001b) for the Trapezium with age of ~ 1.5 Myr and 
IC 348 with - 2.3 Myr (Palla & Stahler 2000), respec- 
tively. From this tendency, we can use the fractions 
of PMS member candidates in order to conjecture the 
age of the W51B clusters. The fractions are 66/168 
(39%), 105/589 (18%), and 45/120 (38%) for G48.9- 
0.3, G49.0-0.3, and G49. 2-0.3, respectively. Subtract- 
ing the field star contaminations, the fractions are given 
by 44/108 (41%), 39/246 (16%), 33/89 (37%) in order. 
These fractions of 20-40% imply that the age of the 
clusters are approximately 2 Myr. It should be the up- 
per limit of the cluster age because we detect only mas- 
sive OB stars whose NIR excess fraction is expected to 
be lower than the low-mass stars due to more rapid disk 
dispersal times (Lada et al. 2000; Haisch ct al. 2001b). 
However, L-band observation is needed for robust es- 
timation of cluster age because if -band involves risk 
under the influence of nebulous environments (Haisch 
et al. 2000). 

VI. SUMMARY 

We have presented results from deep J-, H-, and 
-Kg-band observations of two fields including three com- 
pact H II regions along W51B GMC. We set constraints 
to select the possible cluster members associated with 
each H II region. We investigate the IMF using J-band 
luminosity and estimate total stellar mass from a well- 
defined comparison IMF through the whole mass range 
(from O.OIMq to the mass of the most massive stars). 
In particular, we find the followings: 

- We compute the reddening slope in the direction 
of W51B in the 2MASS system, by a least-squares 
fit with 3-cr clipping. The resultant color excess 
ratio is Ej- H /E H _ Ks - 2.07. 

- We define the cluster member candidates as the 
sources within the distinct bump in the histogram 
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along J—K s color and the sources showing NIR 
color excess emission in the color-color diagram. 

- The mean visual extinctions for the cluster mem- 
ber candidates arc 17 mag for both G48.9-0.3 and 
G49.0-0.3, and 23 mag for G49. 2-0.3. 

- Each IMF based on J-band luminosity corre- 
sponds to the IMF of a young open cluster NGC 
2264. A universal IMF assumption provides total 
stellar mass of ~ 1.4 x 1O 4 M , which is nearly 
3 times of the observed stellar mass. The overall 
SFEs are 7 % for the cloud including G48.9-0.3 
and G49.0-0.3, and 17 % for G49. 2-0.3, which arc 
much higher than the Galactic average value of a 
few percent. 

- We estimate the dynamic age of each H II region 
to be the order of or younger than 1 Myr, in terms 
of the expansion of Stromgren sphere of a compact 
H II region. The JHK S excess fractions of 20-40% 
imply the cluster age of ~ 2 Myr, which should be 
the upper limit of age in consideration of more 
rapid disk dispersal times for massive stars. 
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